Spin-charge separation in cold Fermi-gases: a real time analysis 
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Using the adaptive time-dependent density-matrix renormalization group method for the ID Hub- 
bard model, the splitting of local perturbations into separate wave packets carrying charge and spin 
is observed in real-time. We show the robustness of this separation beyond the low-energy Luttinger 
liquid theory by studying the time-evolution of single particle excitations and density wave packets. 
A striking signature of spin-charge separation is found in ID cold Fermi gases in a harmonic trap 
at the boundary between liquid and Mott-insulating phases. We give quantitative estimates for an 
experimental observation of spin-charge separation in an array of atomic wires. 



One-dimensional (ID) quantum many-body systems 
have been at the center of theoretical and experimental 
interest for the last two decades. Following the semi- 
nal work of Haldane 0, it has been understood that - 
independent of their bosonic or fermionic nature - the 
low-energy behaviour of ID quantum liquids is univer- 
sally described by the so-called Luttinger liquid (LL) pic- 
ture 0, ■ Probably the most remarkable prediction is 
the phenomenon of spin-charge separation in the case of 
Fermions, i.e. the fact that - at low energy - the excita- 
tions of charge and spin completely decouple and prop- 
agate with different velocities. A definite signature of 
spin-charge separation requires the observation of the two 
corresponding branches of excitations in the single parti- 
cle spectral function 0|. In condensed matter systems, 
numerous experiments have looked for spin-charge sepa- 
ration e.g. via photoemission from ID metallic wires on 
surfaces ]J, in ID organic wires in carbon-nanotubes 
0, and in quantum wires in semiconductors, where the 
singular nature of the spectral functions associated with 
spin-charge-separation was observed in tunneling exper- 
iments [7|. In the last few years, ultracold gases in opti- 
cal lattices are providing an entirely new area of physics 
where strong correlations can be studied with unprece- 
dented control and tunability of the parameters. In par- 
ticular these systems open the possibility to investigate 
the transition between three, quasi-two and quasi-one di- 
mension. Recently, an 'atomic quantum wire' configura- 
tion in an array of thousands of parallel atom waveguides 
was realized in ultracold Fermi gases by the application 
of a strong two dimensional optical lattice || . The pos- 
sibility to use cold atoms for studying the phenomenon 
of spin-charge separation was first suggested by Recati 
et al. 0. Their analysis is essentially based on the hy- 
drodynamic Hamiltonian of the LL; the inhomogeneity 
due to the presence of a harmonic trap is treated within 
a local density approximation (see also |loj). In prac- 
tice, with typically less than 100 atoms per atomic wire 
0, observable effects require to use stronger and more 
localized perturbations, where a LL description is not 
applicable. In addition, the effect of boundaries, where 
the local density approximation breaks down, are of cru- 



cial importance. For a quantitative description of spin- 
charge separation in ID cold Fermi gases, it is thus nec- 
essary to use a microscopic description like the Hubbard 
model and properly treat the inhomogeneous case with 
realistic system sizes. Due to the short range nature of 
the interactions between cold atoms, the Hubbard model 
is indeed a perfect description of a situation, in which 
there is an additional optical lattice along the weakly 
confined axial direction (for bosons, the corresponding 
setup has already been realized, see |ll],[l2j). It is an es- 
sential new feature of cold atoms in optical lattices that 
parameters can be changed dynamically and the resulting 
time evolution can be studied. This gives direct access 
to the real-time dynamics of strongly correlated systems, 
a subject hardly studied so far. In this context, a suc- 
cessful method is the recently developed adaptive time- 
dependent density matrix renormalization group (adap- 
tive t-DMRG) 13] which is an efficient implementation 
of Vidal's TEBD algorithm Q in the DMRG framework 
(lif . It has previously been applied to study density per- 
turbations in bosonic ID condensates over a large range 
of interaction strengths |16( . The real-time dynamics of 
interacting spinful Fermionic systems is much harder and 
has not been studied except for very small systems [l7j . 
In our present work, we present numerical results of the 
real-time dynamics of a ID Hubbard model for realistic 
sizes of up to 128 sites. Our main results are: 

(i) real-time calculation showing spin-charge separation 
explicitly in systems of experimentally accessible size 

(ii) the demonstration that spin-charge separation sur- 
vives far outside the low-energy LL regime 

(iii) a quantitative calculation for the effect of spin-charge 
separation at the boundary between a liquid and a Mott- 
insulating (MI) phase which allows to observe the phe- 
nomenon in cold gases without the problems arising from 
the different densities in an array of parallel atomic wires 
and to distinguish experimentally between a Mott- and 
a band insulator. 

Our starting point is the standard Hubbard model 
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for Fermions in ID. Its parameters are the hopping ma- 
trix element J, the on-site repulsion U > between 
Fermions of opposite spin a =1,1 at sites j = 1, . . . , L 
and a spin-dependent local on-site energy £j l(T , describing 
both a possible smooth harmonic confinement and time- 
dependent local potentials. One introduces a 'charge' den- 
sity n c = n^+ni and a 'spin' density n s = nj — njj in a re- 
alization with cold gases, the spin degrees of freedom are 
represented by two different hyperfine levels, and 'charge' 
density is part icle density. Similar to bosons in an opti- 
cal lattice 19], the ratio u = U/J can easily be changed 
experimentally by varying the depth Vq of the optical lat- 
tice. We use units where both J and h are equal to one; 
thus time is measured in units of h/J. In the numerical 
calculations below, we study the dynamics of the Hub- 
bard model using different initial density perturbations 
and the excitations resulting from adding a single par- 
ticle at a given lattice site, which is expected to display 
the same physics as contained in single particle spectral 
functions. Experimentally, the density perturbations may 
be generated by a blue- or red-detuned laser beam tightly 
focused perpendicular to an array of atomic wires, which 
generates locally repulsive or attractive potentials for the 
atoms in the wires. In practice, the perturbations due to 
an external laser field are quite strong, typically of the or- 
der of the recoil energy E r and thus clearly require a non- 
perturbative treatment. In our calculations the length of 
the chains was chosen up to L = 128 sites, keeping of 
the order of several hundred DMRG states. DMRG er- 
ror analysis reveals that all density distributions shown 
here are exact for all practical purposes, with controlled 
errors of less than O(10 -3 ) HJ. 

We start with a homogeneous system which is per- 
turbed by a potential Sj.i localized at the chain center 
which couples only to the f-Fermions, i.e. 

e jA {t) oc exp{-[j — (L — l)/2] 2 /8}r9(-;) (2) 

The potential is assumed to have been switched on slowly 
enough for equilibration, and is then switched off sud- 
denly at time t — 0. In Fig. ^ (a) the density dis- 
tribution of the state at t = 0.2 is shown. The external 
potential J2J) generates a dominant perturbation in the f- 
Fermion distribution by direct coupling and, indirectly, a 
smaller perturbation in the j-density due to the repulsive 
interaction between the different spin species. The wave 
packets in f and j-density hence perform a complicated 
time evolution (Fig. In contrast, the perturbations 
in the spin and charge density split into two wave pack- 
ets each moving outwards. Their respective velocities are 
different as indicated by the arrows in Fig. d(b), sepa- 
rating spin and charge. 

In the limit of an infinitesimal perturbation much 
broader than the average interparticle spacing, both spin 
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FIG. 1: Snapshots of the evolution of the density distribution 
are shown at different times. At t = 0.2, a wave packet is 
present in the center of the system in both the spin and the 
charge density. Each of these splits up into two packets which 
move with the same velocity in opposite directions. The ve- 
locity of the charge wave and the spin wave are different. 



and charge velocities are known analytically from the 
Bethe ansatz |2lj |. To compare our numerical findings 
to the exact charge velocity, we create pure charge den- 
sity perturbations, by applying the potential of Eq. © 
to both species, i.e. e^j = s^ii ancl calculate their time- 
evolution after switching off the potential. The charge 
velocity is determined from the propagation of the max- 
imum (minimum) of the charge density perturbation for 
bright (amplitude n c > 0) and grey (n c < 0) perturba- 
tions, respectively. In Fig. |21 the charge velocities for 
various background densities no and perturbation ampli- 
tudes n c are shown. We find good agreement, if we plot 
the charge velocity versus the charge density at the max- 
imum (minimum), i.e. n c = rio + rj c . This stays true even 
for strong perturbations r\ c ~ ±0.1 which corresponds to 
20% of the charge density. The charge velocity is thus ro- 
bust against separate changes of the background density 
no and the height of the perturbation n c . By contrast, 
the velocity of a spin perturbation varies strongly with 
its height. A possible reason for this may be the nonlin- 
ear cos (^-contribution in the LL description of the spin 
density field, which is only marginally irrelevant, giving 
rise to a nonanalytic contribution to the spin-density re- 
sponse |3j. Nevertheless, our numerical results show that 
for decreasing height of the spin perturbation, its veloc- 
ity approaches the value for the spin velocity obtained by 
the Bethe ansatz. 

In order to compare the behaviour of the density per- 
turbations with that of a single particle excitation, the 
time evolution of the system with one additional particle 
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FIG. 2: Exact results for the charge velocity obtained by the 
Bethe ansatz are (lines) compared to the numerical results of 
the adaptive t-DMRG. The numerical results correspond to 
different heights of the perturbations at various charge back- 
ground densities no. n c is the charge density at the maxi- 
mum/minimum of the charge density perturbation. The un- 
certainties are of the order of the size of the symbols and stem 
mainly from the determination of the velocity. 
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FIG. 4: Time-evolution of charge and spin density pertur- 
bations in the presence of a parabolic trapping potential 
e a ,j = —Vta 2 (j — L/2 + 0.5) 2 E r . MI marks the approximate 
MI region in the absence of the perturbation. The line de- 
noted by idealized is a sketch of the charge density distribu- 
tion without the perturbation. The presence of the perturba- 
tion enlarges the region in which the charge density is locked 
to n c = 1. The arrows show the approximate place of the 
spin perturbation, and the shaded region marks the region 
over which the densities are averaged (cf. Fig. 0. 



FIG. 3: Snapshot of the time-evolution of the charge and spin 
densities of a single particle excitation created at time t = 
at site j = 37 is shown for t = 7.2. 



added at time t = on site j to the ground state, is cal- 
culated numerically. In Fig. [3] a snapshot of the resulting 
evolution of the densities is shown for time t — 7.2. Re- 
markably, as in the case of the density perturbation sep- 
arate wave packets in spin and charge can be seen. This 
demonstrates the phenomenon of spin-charge separation 
directly in a single particle excitation, in close analogy 
to the situation of an inverse photo-emission experiment 
for the single particle spectral functions 0] 

In a specific experiment with arrays of parallel atomic 
wires, it is necessary to take into account that there is an 
additional harmonic trapping potential. Moreover, indi- 
vidual wires have slightly different fillings, which leads to 
an inhomogeneous broadening due to the resulting differ- 
ence in velocities. In order to observe an unambiguous 
signal of spin-charge separation in such a situation, we 
suggest an experimental setup, which relies on the co- 
existence of a MI state and a liquid state in spatially 
separated regions of the parabolically confined system 
|22j . The idea is to use the very different behaviour of 
the charge and spin degrees of freedom in the MI phase. 
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FIG. 5: Time evolution of the charge and spin density 
summed over the sites jo = 25 to ji = 35. Vertical lines 
correspond to the times of the snapshots in Fig. 0] 



In this phase the charge excitation spectrum has a gap, 
whereas the spin dispersion is still linear for small mo- 
menta, and the spin velocity is finite. By contrast, in the 
liquid phase both excitation spectra are linear for small 
momenta. To exploit this, assume the system of one- 
dimensional wires is prepared in such a way, that a MI 
region is present in the center, where the charge density 
is locked at half- filling, n c = 1. At the boundary of this 
MI region liquid regions appear. A localized potential in 
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the liquid region will then create spin and charge den- 
sity waves. Calculated snapshots of the time evolution 
in such a situation are shown in Fig. Evidently, the 
spin density wave propagates into the MI region whereas 
the charge density perturbation is almost completely re- 
flected due to the charge gap in the MI. The presence 
of spin density oscillations which are due to the antifer- 
romagnetic coupling induced by the interaction obscures 
the exact evolution of the spin perturbation. However by 
averaging over several lattice sites - as is always necessary 
in an experiment - the effect of spin-charge separation is 
clearly visible. In Fig. [S]examples for the evolution of the 
sum of the charge and the spin number of particles be- 
tween site 25 and 35, N c and N s respectively, are shown. 
It is clearly seen that the sum of the charge occupation 
does not change, whereas the spin occupation shows the 
moving wave packet. The average spin velocity can be 
determined from Fig. [5] if the distance between the lo- 
calized potential which generates the perturbation and 
the center of the region over which the density is mea- 
sured is known. Here, the spin velocity is found to be 
v s w 1.1 J /% which agrees nicely within the expected ac- 
curacy with the value of v s (n c = 1) = 1.2J/H of the 
Bethe ansatz. The very different propagation behaviour 
of charge and spin can as well be used experimentally 
to distinguish between a MI and a band insulator: In a 
band insulator not only the velocity of the charge, but as 
well of the spin would vanish, whereas, as used above, in 
the MI the spin velocity stays finite. 

In order to quantify the requirements for an experimen- 
tal observation of spin-charge separation in cold Fermi 
gases, we finally discuss typical parameters which need 
to be achieved in a setup with an array of atomic wires 
0- Such an array consists of several thousand parallel 
wires with typically less than 100 40 K atoms each. In 
addition to the smooth axial confinement potential with 
frequency to z w Ef/N rj 2ir ■ 275Hz (corresponding to 
Vf w —0.0035), realization of a ID Hubbard model re- 
quires adding a strong periodic potential along the tubes. 
For 40 K and a standard lattice constant a = 413nm the 
recoil energy is E r « 7kHz. An optical lattice of height 
Vo = 15E r then gives an on-site repulsion U ~ 0.17 E r , 
where we have used a standard value for the s-wave scat- 
tering length a s ~ 174a n fo r the F = 9/2 rrif = —9/2 
and mf = —7/2 states 23]. The resulting dimension- 
less interaction u ~ 22 then leads to a central MI region 
with a typical size of around 20 sites. With this param- 
eters, the time in which the spin wave travels 20 sites 
is of the order of a few ms. The creation of state se- 
lective potentials for two different hyperfine states may 
be done by using laser light whose frequency falls be- 
tween the respective transitions. This might be difficult 
for the F = 9/2 m f = -9/2 and m f = -7/2 levels, but 
should be possible - for instance- using the F = 9/2 and 
F = 7/2 levels. The l/e 2 -radius of the potential (Eq. [2J 
is taken to be four lattice sites, which could be realized 



approximately by a laser of an l/e 2 -radius of 2.1/jm or 
less. Finally, to ensure that finite temperature does not 
destroy the Mott insulating behaviour by thermal acti- 
vation, the energy scale ksT should be smaller than the 
Mott energy gap. Already the very first experiment of 
ID fermions in an optical lattice is very close to matching 
those conditions. Improvements can be reached reducing 
the axial confinement frequency. 

To conclude we have performed numerical simulations 
of the time-evolution of charge and spin density pertur- 
bations in the ID Hubbard model. We clearly observe 
the separation of spin and charge as a generic feature of 
ID fermions, far beyond the low-energy regime where a 
Luttingcr liquid description applies. In addition, an ex- 
periment is suggested which exhibits the separation of 
the two modes from the perfect reflection of density ex- 
citations at the boundary to a Mott insulating state 

We gratefully acknowledge very fruitful discussions 
with Th. Giamarchi and with the experimental group 
of T. Esslinger, in particular H. Moritz. 
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